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Abstract 

We investigate the dynamics of impurities walking along a semi- 
conductor layer assisted by thermal noise of strength D and external 
harmonic potential V{x). Applying a nonhomogeneous hot tempera- 
ture in the vicinity of the potential minimum may modify the external 
potential into a bistable effective potential. We propose the ways of 
mobilizing and eradicating the unwanted impurities along the semi- 
conductor layer. Furthermore, the thermally activated rate of hopping 
for the impurities as a function of the model parameters is studied in 
high barrier limit. Via two state approximation, we also study the 
stochastic resonance (SR) of the impurities dynamics where the same 
noise source that induces the dynamics also induces the transition 
from mono-stable to bistable state which leads to SR in the presence 
of time varying field. 
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1 Introduction 



Recently the physics of semiconductors has received considerable attentions 
as devices made from semiconductors are vital in the construction of modern 
electronic equipment. One of the beneficial features of these semiconduc- 
tors is the possibility to adjust their conductivities by introducing impurities 
(dopants) into their crystal lattice. The level of conductivity can be con- 
trolled by the amount and type of impurities. Unlike conductors, their con- 
ductivity increase with temperature; this property makes them even more 
important. 

Since the conductivity of the semiconductor relies not only on the con- 
centration of the dopants but also on the thermal background strength of the 
medium, the physics of thermal diffusion of the impurities along the semi- 
conductor layer has attracted considerable attentions [11 |2l [3l H] . Recently, 
the proposal by [1] unveiled the method of eliminating unwanted impuri- 
ties from the region of the semiconductor via movable external harmonic 
potential along the semiconductor. The same group has investigated how 
to control the diffusion of impurities via an external potential which has an 
advantage of preserving the crystal structure as diffusion can take place at 
low temperature [2] . Their numerical results exhibit that the diffusion of the 
dopant increases when the strength of the external potential decreases. At 
high doping density, the internal field becomes significant and it renormal- 
izes the effect of the external potential. Thus the diffusivity of the dopant 
increases as the dopant density increases. 

Because manipulating the diffusion of impurities (either acceptors or 
donors) to a desired region of the semiconductor layer is vital, in this work we 
propose different ways of mobilizing the impurities along the semiconductor. 
We consider the impurities walking from one lattice trap to the other lattice 
trap with a trap depth potential along the semiconductor layer assisted by 
thermal noise of strength D and the external harmonic potential. Further- 
more, a nonhomogeneous hot temperature is applied in the vicinity of the 
potential minimum which may modify the external potential into a bistable 
effective potential. Neglecting the interaction between the dopants at low 
impurity density, we show that the diffusion of impurities can be controlled 
not only by varying the external potential but also by altering the intensity 
of the temperature of the hot spot and trap depth of the impurities. We 
also propose the way of eradicating the unwanted impurities from the cer- 
tain region of the semiconductor by varying the different model parameters 
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without considering movable external potential. Furthermore, we investigate 
the dependence of the rate on different model parameters. 

One very crucial but unexplored issue is the way of enhancing the mo- 
bility of the impurities along the semiconductor. This can be achieved by 
applying symmetry breaking fields such as time varying signal which may 
lead the system into stochastic resonance. Stochastic resonance has been 
widely studied over the last few decades [SI [3 El El [ID] and has been adopted 
to describe an interesting statistical property of periodically modulated and 
noise driven multi-stable dynamical system. It exhibits that under proper 
condition an increase in the input noise level results in an increase in the 
output signal-to- noise ratio. 

The phenomenon of SR was first introduced by Benzi [S]. Later, the 
idea of stochastic resonance has been implemented in many model systems. 
Notable examples include: SR for confined systems |Tl], SR for complex 
systems such as polymers [71 El [121 [IS] and SR for a Brownian particle moving 
across a porous membrane [T^. The vast majority of studies on stochastic 
resonance have been focused on analyzing the dynamics of a bistable system. 
Among such models, the two state model [51 [IS] has been proven to be 
extremely useful in the understanding of the stochastic phenomenon offering 
a simple framework able to provide analytical results 

In this paper, via two state approximation, we study the stochastic res- 
onance (SR) of a particular noise sustained dynamics where the same noise 
source that induces the dynamics also induces the transition from mono- 
stable to bistable state which leads to SR in the presence of time varying 
field. We find that the spectral amplification rj attains an optimal value at 
a certain finite value of noise strength Dopt- As the temperature of the hot 
spot increases, the peak of rj decreases. On the other hand, the value of 
thermal strength Dopt increases as the intensity of the nonhomogenous hot 
temperature increases. The strength of trap depth also considerably affects 
the rj. We find that the peak of the rj decreases as the strength of the poten- 
tial trap increases. Furthermore, we explore the dependence of rj on other 
model parameters. 

It is important to note that the results and model presented in this work 
are not specific to the impurities dynamics. Rather the present study serves 
as a basic paradigm in which to understand diffusion and noise induced 
nonequilibrum phase transition in discrete systems. Thus the present model 
is of broad interest in various fields. 

The rest of the paper is organized as follows: in Section 2, we present 
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the model. In Section 3, we study the diffusion of the impurities in the 
presence of a nonhomogeneous temperature. In section 4, we study the SR 
of impurities and explore how the t] for the impurities behaves as a function 
of the model parameters. Section 5 deals with summary and conclusion. 



2 The model 

We consider a one-dimensional system where non-interacting impurities jump 
from one lattice trap to the other lattice trap assisted by thermal noise and 
external potential energy which is monostable 

V{x) = Vox'' (1) 

where Vq and x denote the potential energy and the position of the impurities, 
respectively. When background temperature is homogeneous, particles con- 
centrate around the potential minimum. On the other hand, in the presence 
of a nonhomogeneous temperature background, 

T{x)=T, + nexp[-^] (2) 

which is hot around the potential minimum and decays to a lower tempera- 
ture as one goes away on both sides from the potential minimum, the system 
may undergo a phase transition in which the particles pile up around two 
points of the potential minima. Here the parameters and designate the 
temperature of the hot and cold reservoirs, respectively, while a denotes the 
standard deviation. 

At low impurity density, the impurity dynamics is governed by [171 [18] 

dP{x,t) d , V'{x) ^ ^ ,r.^ ^ d ^ * ^, 

at ox kBT[x) kBT[x) ox kBT[x) 

(3) 

where P{x, t) is the impurity density at a position x and time t. The im- 
purities jump from one lattice trap, which has an internal potential depth 
(trap depth) $, to the next lattice trap along the one dimensional lattice. 
The steady state probability distribution of the impurity Pss{x) is given by 

Pss{x) = CeT^expl- / TT^dy] (4) 

Jo kBT[y) 
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where C is a normalizing constant. Rearranging Eq. 4, one gets the proba- 
bihty distribution 

Pssix) = Coe (5) 



with an effective potential energy Veff{x) which is given by 



+ a 



[1 + ae 



x2 



(6) 



where dimensionless parameter a relates Th and Tc as Th = aTc. The effective 
potential energy is bistable when the condition 



$ > 2Voa' 



a + 1 



a 



(7) 



is met. In this regime, the location of the saddle point is at x = while the 
two symmetric stable points are located at 



Xr, 
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The effective potential energies at the saddle and stable points are given as 



VeffiO) 



l + a) 
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and 



-$ + 2ctVo 



1 + ln 
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+ In 



[1 + a){^ - 2a^Vo 



(10) 

respectively. The two minima of the potential are separated by a barrier of 
height AV^-^^ = Ve//(0) — Ve//(a^m)- The curvatures at the barrier top Uq 
and the well minima take simple forms 



[-a^/a^ - 2a{l + a)Vo + 2(1 + afVo) 
(1 + aY 



and 



|^($-2aVo)^ln[a 



2a'^Vn 



1)], 



(11) 



(12) 



respectively. 

Transforming Eq. (4) into Eq. (5) is equivalent to converting the non- 
equilibrium problem to an equilibrium one where the usual Boltzmanns 
statistics holds. Hereafter we consider impurities undergoing a random walk 
motion along the effective potential assisted by the thermal kicks D = ksTc 
where is the Boltzmann constant. One can note that the information 
about the non-equilibrium features of the system is stored in the effective 
potential. In the next section we explore how the dynamics of the impuri- 
ties under the effective potential V^-^-^ behaves as a function of the model 
parameters. 

3 Impurity diffusion in a nonhomogeneous tem- 
perature 

We consider non-interacting impurities (either acceptors or donors [2] ) of low 
density hopping from one lattice site to the other site within the semicon- 
ductor layer. It is assumed that only the donors are sensitive to the external 
potential and hereafter the term impurity refers to the donor only. Expos- 
ing the impurities to the external harmonic potential compel the particles to 
accumulate around the potential minimum. Furthermore, applying a non- 
homogeneous temperature profile in the vicinity of the potential minimum, 
may modify the external potential into a bistable effective potential as long 
as the condition given in Eq. (7) is obeyed. Otherwise the effective potential 
is a monostable potential. 

Figure la gives a plot of the effective potential V^^^^ versus x for two 
different values of a. The effective potential has two stable potential minima 
at diXm- Exploiting Eq. (6), one can see that the barrier height AV^^^ and 
the width of the effective potential 2xm increase as a and the trap depth $ 
increase. On the contrary, as the strength of the external potential Vq and a 
decrease, the barrier height AV^-^-^ and the width 2xm increase. 

We plot the dependence of the probability distribution Pss{ func- 
tion of X in Fig. lb. The figure depicts that the impurities accumulate 
around the two potential wells. As $ increases, the probability of finding 
the impurities around the two stable points increases. When $ increases, 
not only the peak of Pss{x) intensifies but also the distance between the two 
peaks 2xm increases. 
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Figure 1: (a) (Color online) (a) The effective potential V^-^-^l function 
of X (nm) for parameters choice a = 2 (black line), a = 4 (red line), $ = AeV, 
a = l.lnm, = 300K and Vq = OAeV/nrn^. (b) Probability distribution 
Pssix) as a function of x for parameter choice Tc = 300K, a = 1.5, Vq = 
QAeV /nm? and cr = l.lnm. The black and red lines stand for $ = 2eV and 
$ = QeV , respectively. 

Let us now discuss how one can translocate the dopants to a desired lo- 
cation along the semiconductor layer. By looking at how V^^^ and Pss{x) 
behave, when the different model parameters vary, one can infer the po- 
sition of the impurities in the medium. As discussed before when a and 
$ monotonously increase or as Vq and a constantly decrease, the particles 
which are accumulated around the two potential minima leave the central 
region x = and migrate to the peripheral part of the semiconductor layer 
creating a depletion zone around a; = 0. This suggests ways of controlling the 
diffusivity of the particles by tuning these parameters. Further increase in a 
or as one further decrease Vq and a, the impurities may become eliminated 
from a certain region of the semiconductor layer. 

The impurities that are exposed to the effective potential exhibit noise- 
induced nonequilibrium phase transition [15] where the noise plays a coun- 
terintuitive role by inducing ordering phenomena. Exploiting Eq. (7), one 
can see that the regimes of noise induced ordering phase transition are de- 
marcated by the critical values of the different model parameters. The crit- 
ical potential energy Vq and the critical trap depth $* have a simple form: 
V* = ($a/(2a2(a + 1)) a nd $* = 2Voa^ {{a + l)/a) while the other critical 

points are given by a* = y^(a$)/(2Ko(a + 1)) and a* = {2VQa^)/{<^-2VQa^). 
When < $ < $*, a < a*, (7 > (7* or V"o > Vq*, the effective potential is 
a monostable potential and the dopants concentrate around the x = po- 
sition. In such a case the impurity dynamics is similar to the previously 



exposed works [U On the other hand when $ > $*, a > a*, cr < cr* or 
Vq < Vq, the effective potential has two stable minima and the impurities 
split into two impurity-rich regions which coexist with impurity-poor regions. 
In this regime, the donors pile up towards the two potential minima while 
the oppositely charged particles hop in the opposite direction showing that a 
series of p — n—p or n—p — n junctions can be fabricated by manipulating the 
diffusion of the impurities. Note that the effect of internal field is neglected 
which is appropriate at low doping level. 

Now let us study how the thermally activated barrier crossing rate of 
the impurities behaves. Let us consider non-interacting impurities initially 
situated at one of the potential minima. Due to thermal fluctuation, the par- 
ticles cross the potential barrier assisted by the thermal kicks they encounter 
along the reaction coordinate. The crossing rate for the impurities in a high 
barrier limit AV^^^ ^ ksT is approximated 120] as 
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Figure 2: (a) (Color online) Escape rate R versus in a unit of Kelvin 
for parameters choice a = 1.5 (black line), a = 3 (red line), $ = leV and 
a = Inm. (b)The crossing rate R versus in a unit of Kelvin for parameters 
choice Vq = OAeV /nm^ (black line) and Vq = O.SeV /nm'^ (red line). Other 
parameters are fixed as a = 2, $ = 2eV and a = Inm. 

Figure 2a plots the dependence of the rate R on Tc. The particles cross 
the potential barrier at the expense of the thermal background kicks. Thus 
the rate increases with T^. When a steps up, AV'^^'^ and 2xm increase and 
the particles attain difficulty in crossing the high barrier. Hence the rate 
becomes lower with a. Figure 2b exhibits that when Vq rises, the escape 
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rate increases due to the fact that increasing in Vq results in a lower barrier 
height. Exploiting Eq. (13), one can also see that the rate intensifies when 
$ decreases or as a increases. 



4 Stochastic resonance of the impurities 

Let us now study the noise-assisted dynamics for the bistable potential in 
the presence of time varying signal. The interplay between noise and time 
varying force in the bistable system may lead the system into stochastic 
resonance as long as the random kicks are adjusted in an optimal way to the 
recurring external force. Next we study the dependence of the SR on the 
model parameters employing two state approximation. 

Via two state model approach [6], [16], two discrete states x{t) = ±Xm 
are considered. Let us denote n+ and n_ to be the probability to find the 
impurity in the right (x^) and in the left (— x^) sides of the potential wells, 
respectively. In the presence of time varying weak periodic force (AC field) of 
type v(t) = AQCos{flt), the master equation that governs the time evolution 
of n± is given by 

n±it) = -W±it)n± + W^{t)n^ (14) 

where W+{t) and W-{t) correspond to the time dependent transition proba- 
bilities towards the right (x^) and the left (— x^) sides of the potential wells, 
respectively. The time dependent rate [SI [IS] takes a simple form 



W± = Rexp 



: ™" Aq cos (Qt) 



(15) 



where R is the Kramers rate for the particles in the absence of periodic force 
Aq = 0. For sufficiently small amplitude, one finds the spectral amplification, 
r/ to be P [21] 

The input signal that modulates the symmetric bistable system makes 
one stable state less stable than the other, alternatively, over half forcing 
period. When the random switching frequency matches the forcing angular 
frequency by tuning the noise intensity, the system attains the maximum 
probability of escaping out of the less stable state into the more stable one, 
before a random back switching event takes place. Due to this reason, the 
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particle is more likely to be in the more stable state. Consequently rj attains 
a maximum value at a particular noise level {Dopt)- Let us now explore how 
the rj behaves as a function of the trap depth $, standard deviation a, applied 
potential energy Vq and a. 
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Figure 3: (a) (Color online) Spectral amplification t] versus in a unit 
of Kelvin for parameters choice Vq = OAeV /nm^ (black line) and Vq = 
O.SeV /nm'^ (red line). Other parameters are fixed as a = 2, $ = 2eV and 
a = Inm. (b) Spectral amplification rj versus in a unit of Kelvin for 
parameters choice a = 1.5 ( black line), a = 3 (red line), $ = 2eV and 
a = Inm. 

Figure 3a plots the dependence of the rj on the noise strength for 
parameter choice Vq = OAeV/nm'^ for black line and Vq = O.SeV/nm^ for 
red line. Other parameters are fixed as a = 2, $ = 2eV and a = Inm. 
In the limit Tc — and Tc — )■ oo, r/ — )■ 0. This is because when the noise 
intensity Tc is too small, the intrawell crossing rate is too small while for 
large Tc a similar loss of synchronization occurs since the particles flip too 
many times between the two stable points. In between t] attains an optimal 
value at optimal T°^*. The peaks of t] increases with Vq. On the other hand, 
T°^* increases as Vq decreases. Figure 3b exhibits the dependence of t] on Tc 
for parameters choice a = 1.5 for black line, a = 3 for red line, $ = 2eV, 
a = Inm and Vq = O.AeV/nm"^. When a increases, T°p* increases; the peak of 
the rj rises when a diminishes. This loss of synchronization with the increase 
in a is plausible because as a increases, AV^J^ and the width 2x™ mcrease; 
the particle crosses the high potential barrier at the expense of higher thermal 
strength Tc. 

In Fig. 4a, we show how the rj behaves as $ varies. When $ increases, 
the peak of t] becomes smaller. T"^* shifts to the right as $ increases showing 
the response of signal to the background noise strength is significant at lower 
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Figure 4: (a) (Color online) Spectral amplification t] versus Tc in a unit of 
Kelvin for parameters choice $ = 2eV (black line) and $ = 2.5eV (red line). 
Other parameters are fixed as a = 1.5, Vq = OAeV/nm'^ and cr = Inm. (b) 
Spectral amplification rj versus Tc in a unit of Kelvin for parameters choice 
1 (black line), a = 1.1 (red line), $ = 2eV, a = 2 and Vq = OAeV/nm'^. 
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values of $. On the other hand, as a monotonously decreases, the resonance 
becomes more significant. When a diminishes, T°^* increases (see Fig. 4b). 

These results reveal the weak signal passing through the semiconductor 
layer can be amplified and detected by tuning the temperature of the hot 
locality since the noise level can be directly affected by the hot temperature. 
By tuning the angular frequency Q, one can control the mobility of the 
impurities to a desired location. At the resonance temperature T"^*, the 
impurities undergo a fast unidirectional drift from the less stable potential 
minima to the more stable one over half forcing period revealing novel way of 
achieving a fast transportation of impurities along the semiconductor layer 
without exposing the impurities to a higher temperature. Thus the present 
study is crucial in the designing of artificial semiconductor. 

At this point we stress that this theoretical work can be realized exper- 
imentally. As shown in the schematic diagram (Fig. 5), one can model the 
harmonic potential via the metallic gates which are kept at a certain voltage 
[H 12]. The metallic gates are situated at the top of the sample semiconduc- 
tor. Similar to the work [H [2] , the harmonic potential is extended along the 
X direction. A tiny region around the potential minimum x = is heated 
up and this creates a depletion zone of the impurities around x = 0; the 
hot locality forces the impurities to migrate towards the peripheral regions. 
The diffusion of the particles can be controlled by tuning the different model 
parameters. Applying periodic signals such as AC field along the x direc- 
tion, modulates the effective potential. With the proper adjustment of the 
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Figure 5: (Color online) Schematic diagram showing a semiconductor which 
is exposed to an external potential (in x direction) and a nonhomogeneous 
temperature background, which is hot around the potential minimum (x = 0) 
and decays to a lower temperature as one goes away on both sides from the 
potential minimum. The cross section along the z direction is not shown. 
The impurities are assumed to be uniformly distributed along the y and z 
directions. Furthermore, a time varying field (AC field) is applied along the 
X direction. 



different parameters, the system may show the stochastic resonance and this 
arrangement could be used to amplify very weak signals. 

5 Summary and conclusion 

This theoretical work exposes the way of manipulating the disfusiblity of the 
impurities along the semiconductor layer by locally heating the semiconduc- 
tor layer around the potential minimum of the exerted external harmonic po- 
tential. The theoretical results obtained in this work depicts that the dopant 
mobilizes to the peripheral regions when the trap depth $ and a increase or 
when a and the external potential Vq decrease. The thermally activated rate 
for the impurities is also studied at high barrier limit AV^^^ ^ ksT. It is 
shown that the rate increases with Vq and a, and it grows smaller when a 
and $ increase. 
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In the presence of periodic signals, the dependence of the spectral ampli- 
fication rj on the different control parameters is explored. The peak of rj rises 
when $ and a monotonously fall. When a and Vq increase, the response 
of signal to the background temperature is significant. The magnitude of 
optimum noise intensity T°^* becomes considerable with $ and a, and when 
Vq and a grow smaller. 

The model presented in this work is not limited to the impurities dy- 
namics. Rather the present model is of broader interest in various fields and 
serves as a basic paradigm in which to understand diffusion and noise in- 
duced nonequilibrum phase transition in discrete systems. In conclusion, the 
proposed model is crucial in designing artificial semiconductors. We believe 
that a semiconductor device under the suggested model can be designed to 
detect weak signals of extremely small modulation. 
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